Perturbative analysis of topological entanglement entropy from conditional 

independence 



Isaac H. Kim^ 

^Institute of Quantum Information, California Institute of Technology, Pasadena CA 91125, USA 

(Dated; December 6, 2012) 

We use the structure of conditionally independent states to analyze the stabihty of topological 
entanglement entropy. For the ground state of quantum double or Levin-Wen model, we obtain a 
bound on the first order perturbation of topological entanglement entropy in terms of its energy 
gap and subsystem size. The bound decreases superpolynomially with the size of the subsystem, 
provided the energy gap is nonzero. We also study the finite temperature stability of stabilizer 
models, for which we prove a stronger statement than the strong subadditivity of entropy. Using 
this statement and assuming i) finite correlation length ii) small conditional mutual information of 
certain configurations, first order perturbation effect for arbitrary local perturbation can be bounded. 
We discuss the technical obstacles in generalizing these results. 



I. INTRODUCTION 

Topological order is a new kind of order that cannot 
be described by Landau's symmetry breaking paradigm. 
Properties of these exotic phases include a ground state 
degeneracy that depends on the manifold, anyonic statis- 
tics, and long range entanglement .^^K^ Such phases are 
expected to be stable against generic perturbation if its 
strength is sufficiently weak and its interaction range is 
bounded. Indeed, it was shown by several authors that 
the spectral stability follows under a set of reasonable 
assumptions 

If the energy gap remains open under the perturba- 
tion, one can adiabatically continue from the ground 
state of the original hamiltonian to the ground state of 
the perturbed hamiltonian.'^ Since the generator of this 
flow consists of quasi-local terms which decay almost ex- 
ponentially, the perturbed hamiltonia n has similar prop- 
erties to the unperturbed hamiltonian.'^'^'^ For example, 
one can define local operators that create defects with 
well-defined energies and string operators that can move 
around the defects. One may argue that the long range 
entanglement in the ground state can be preserved in 
a similar vein, although one must define precisely what 
long range entanglement is. 

Long range entanglement in a 2D system refers to the 
nontrivial constant subcorrection term of the entangle- 
ment entrop y, also known as the topological entangle- 
ment entropy.SEEfi) While a proof with full mathematical 
rigor has not been established to the best of author's 
knowledge, it is widely accepted by now that topological 
entanglement entropy is a universal constant that char- 
acterizes the phase of the gapped quantum many-body 
system. If one accepts the topological quantum field the- 
ory description of the low energy physics, there is a sim- 
ple explanation as to why the topological entanglement 
entropy remains stable against generic perturbation.Sl 
There are also mounting numerical evidences suggesting 
the stability.'iStll! 

Presence of the long range entanglement can be in- 
terpreted as a consequence of some nontrivial nonlocal 



constraint. For example, in the ground state of a 2D 
gapped system supporting anyonic quasiparticles, total 
charge enclosed in some region must add up to be a trivial 
charge. However, the existence of the constant subcor- 
rection term alone does not necessarily imply that the 
nature of the constraint is quantum. 3D toric code at 
finite temperature has nonlocal contributions to the en- 
tanglement entropjC^, yet such state can be mapped to a 
Gibbs state of a classical hamiltonian under local unitary 
transformation.^^ We wish to understand if this nonlocal 
contribution to the entanglement entropy is an invari- 
ant of the phase. We would also like to understand the 
mechanism behind their stability, instead of arguing on 
the ground of effective field theory. In such pursuit, we 
introduce a property of these states that has apparently 
been unnoticed so far with few notable exceptions. 

The key property is the conditional independence. Tri- 
partite state pabc is referred to be conditionally in- 
dependent if its conditional mutual information I [A : 
C\B) = Sab + Sbc — Sb — Sabc is 0. Hastings and 
Poulin showed that even for a system with long range 
entanglement, there exists a set of subsystems that are 
conditionally independent.'^ To see this, note that the 
entanglement entropy of a topologically ordered system 
can be expressed as Sa = a\dA\ — ja, where \dA\ is 
the boundary area and ja is a constant that only de- 
pends on the topology of A. If the topologies of AB, 
EC, B, and ABC are all identical, 7 as well as the 
boundary contributions cancel out each other. Proving 
such statement for a generic quantum many-body system 
is a hard problem. However, the entanglement entropy 
formula for the g roun d state of some exactly solvable 
models are known.l^ESl Por such systems, the conditional 
independence can be easily shown. Conditional indepen- 
dence is also the key ide a beh ind quantum belief prop- 
agation(QBP) algorithm .'i^Hi Success of the QBP algo- 
rithm indicates that the conditional mutual information 
for certain configuration is like ly to be small for noncrit- 
ical finite temperature systems.' ^^ * ^° ' 

Main message of the present paper is that the condi- 
tional independence of certain subsystems strongly con- 
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strains the structure of these states so as to ensure the 
robustness of the nonlocal quantum correlation. Admit- 
tedly our result is limited to either i) the first order of the 
perturbation series of the exactly solvable models or ii) 
the perturbation that has a special structure. However, 
generalizations to both higher orders of perturbation se- 
ries and more general models seem to be hindered by our 
lack of understanding about approximately conditionally 
independent states rather than that of the physical prop- 
erties of such systems. 

It has been known in quantum information community 
that the structure of states that are conditionally inde- 
pendent is significantly constrained compared to that of 
the the general quantum state.l^^lt^ll particular, exact 
conditional independence implies that subsystems form 
a quantum Markov chain. This structure, together with 
the locality of the hamiltonian, reveals why topological 
entanglement entropy changes very little, at least up to 
the first order of the perturbation series. A statement 
that extends to the approximate conditional independent 
states are not known. In fact, a classical statement that 
relates conditional mutual information to a Markov chain 
is known to be false for quantum states.f^^ 

Rest of the paper is structured in the following way. 
In Section ITTl we introduce the notations. Section ITTT 



explains the technical tools. We apply these tools in Sec- 
tion |IV| and |V] which are respectively dedicated to the 
zero temperature and the finite temperature states. We 
discuss technical obstacles in generalizing the results to 
both higher orders and more general models in Section 
\VJ\ We discuss the implication of these results and some 



open problems in Section VII 



This approximation is motivated from the fact that a 
correlation generated by local hamiltonian falls off expo- 
nentially outside an effective lightcone. The quasilocal 
operators generated by such time evolution can be ap- 
proximated by a local operator supported on a ball of 
finite radius R, with the correction term decreasing su- 
perpolynomially with R. 




^ : correlation length 



FIG. 1: The shaded region represents an effect of the pertur- 
bation that is smeared out in space. We shall approximate 
this effect by a strictly local operator with a finite radius R. 
The correction decreases superpolynomially with R. 



Entanglement spectrum of a subsystem A is defined as 
Ha = —Ia^ ^ log PA, where pA is the reduced density 
matrix of A. We define conditional mutual spectrum as 



H 



A:C\B 



H 



AB 



H 



BC 



Hn-H 



ABC- 



Note that 



Ti{pabcHa:C\b) = I{A : C\B) 



(2) 



II. NOTATIONS 

The Hilbert space has a tensor product structure ®iHi 
where 'Hi corresponds to the local Hilbert space located 
at vertices of a square lattice. Local Hilbert space dimen- 
sion is d. We assume periodic boundary condition with 
sufficiently large system size. We define a set of operators 
having nontrivial support on T-La as B{'Ha)- The Bound- 
ary of subsystem A is denoted as dA. \A\ represents the 
volume of A and similarly \dA\ is the boundary area of 
A. We set the size of the subsystems to be 0{l) unless 
specified otherwise. 

We consider a family of hamiltonian H{s) — Hq + sV 
and study its behavior in the vicinity of s = 0. Both the 
original hamiltonian Hq = hi and the perturbation 
V = J2i ■^i consists of sum of terms that are supported 
on a ball of radius tq and the interaction strength is uni- 
formly bounded by J, i.e. \\hi\\, \\vi\\ < J. || • • • || is ^oo 
norm. We denote the spectral gap as r(s). 

Following Bravyi et al.'s constructiorP^, we define an 
approximation of a quasi-local operator as follows. 

[0]a = ^r^Jr a40)^Ia^ (!) 



We also define (• • • ) — Tr(p • • • ) as an expectation value. 
Throughout the paper, constants c and c' denote numer- 
ical constants, and their exact values may be different in 
each contexts. 



III. TECHNICAL TOOLS 

We have extensively used three technical tools in the 
present paper. The ideas that motivate each of these 
tools are as follows. First, local perturbati on p erturbs 
locally with superpolynomially decaying tail.l^ Second, 
perturbation that acts locally can be bounded by us- 
ing deformation moves. Deformation move refers to a 
chain rule of conditional mutual spectrum, analogous 
to the chain rule of conditional mutual information. 
Third, effect of the superpolynomially decaying tail can 
be bounded by regularizing the entanglement spectrum. 

The locality estimates come from Lieb-Robinson 
bound. The deformation moves and the regularization 
of the entanglement spectrum is a more refined treatment 
of the tools introduced in Refi28l 
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A. Regular izat ion of entanglement spectrum 

We extend some of the results in RefHH) 

Definition 1. Regularized entanglement spectrum 
with a cutoff A is 



Hi 



p>l/A 



(3) 



where \i) is an eigenstate of pA with an eigenvalue pi 
Lemma 1. 

lMBA^|i<^logA, 



where A^^ Ha~ H^- 

Proof Purify pab to |'0)abc- 
ing Schmidt decomposition. 



'ABC 



admits the foUow- 



'ABC 



i=l 



A IVBC ' 



(5) 



where piS are eigenvalues of pA and (N)sc) 
basis states for the Hilbert space HAiT^Bc)- 

For any operator O E B{T-Lab)i it allows the following 
decomposition. 



^ = E E d-^Tr(C/A,. Ub,,0)uI, ® Ul^, (6) 



where UA,i{UB,j) are unitary operators that are sup- 
ported on A{B) with appropriate normalization condi- 
tions. 



Tr{UA.,ulj) = dAS,^ 



(7) 



In other words, {UA,i/^/dA} {{U b ,i / ^/dn}) is a com- 
plete set of orthonormal basis for B{'Ha) {B{T-Lb)) un- 
der Hilbert-Schmidt inner product (Oi, O2) = Tr(0|02)- 
Such basis set always exists for a finite dimensional 
Hilbert space.l^^l Equation [g] is equivalent to the following 
expression. 



o = J20b,^uI^ 



i=l 



where 



Os,,: = —TrA{UA,^0) 
dA 



t.d^^'(UA,^UB,0)Ul^. 



(8) 

(9) 
(10) 



Also, OB,i can be bounded as follows. 

-T- sup X(0lBc(«lAt^^:»O|i; 
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_B.il 



dA 



dA 



^E;rll^-^'''^ll = ll^l 



a\4>)bc (11) 

(12) 



Rewriting TY(pABA'XOB,i U\ J 



I ABC^ 



(VUsc ^AOB,^ ® f/I,, |V')Ai3c = TrlplA^C/^ .plog,,), 

(13) 



(4) where OB,i — VOB.iV^ with isometry V — J2i N 



O"^ is the transpose of O. Equation 13 can be bounded 

by 



l4AlM|{//piof||<^logA||0, 
A2 



Summing over all i, we get 



|Tr(p^sA^O)|<||0||-4logA 
As 



(14) 

(15) 
□ 



Following corollaries can be easily proven by a judi- 
cious choice of A. 

Corollary 1. 

\Tr{pAB log paO)\ < 6||0||logdA (16) 

Corollary 2. Consider a connected correlation func- 
tionC{Oi,02) = (O1O2) - (Oi)(02). //C(0i,02) < 
||Oi|l||02|ie/or all Oi,02, 

|C(7?^,0)| <e||0||(181ogd^-^41og-). (17) 

e 

B. Deformation moves 

Author has introduced a set of deformation moves to 
show that C{Ha:C\b,0) can be bounded for an operator 
O supported on one of the subsystems, provided that i) 
area law holds approximately ii) correlation decays suf- 
ficiently fast iii) certain information-theoretic conjecture 
is correct 

Here we construct a similar, yet slightly different state- 
ment. As in Ref,28, the statement concerns a correlation 
bound between Ha:c\b ^-ud an arbitrary operator O. The 
difference is that here we relax the condition on the sup- 
port of O: O is allowed to be located anywhere, as long 
as its support is sufficiently small compared to the sub- 
system. 

The price we have to pay is the following. Instead 
of relying on an information-theoretic conjecture that is 
expected to hold for any quantum states, we impose a 
condition on the reduced density matrices. 
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Definition 2. pabc is CQ-bounded if 



\Trc{p 



abcHa:c\b, 



\i<coIiA:C\B). 



(18) 



Note that all classical states are 1-bounded. Reduced 
density matrices of finite temperature Gibbs state for the 
so called "stabilizer models" are also 1-bounded. De- 
tailed explanation about these states shall be presented 
in Section|vj If I{A : C\B) = 0, conditional mutual spec- 
trum is 1-bounded by Petz's theorem. More specifically, 
Petz showed that 



Ha:C\1 







(19) 



if and only if I {A : C|B) = O.SH 

Following Ref l28L given a conditional mutual spectrum 
Ha:C\Bj we shall refer B as a reference party. A and C 
shall be referred as target parties. Diagrammatically the 
reference party will be denoted with a 'R' sign and the 
target parties will be denoted with 'T' signs. 

The key idea behind the deformation move is that for 
any local operator O, one can decompose Ha:c\b into 
HAi-.CilBi such that either i) I{Ai : Ci\Bi) = o(l) or ii) O 
is sufficiently far away from AiBiCi. Such decomposition 
can be expressed as a linear combination of the following 
chain rule, which can be verified easily. 



H 



AiA2:C\B 



= H 



A2:C\B 



H 



Ax:C\A2B, 



(20) 



Nevertheless, we found it instructive to define three ele- 
mentary deformation moves to explain this technique. 

The first step in the deformation procedure is to apply 
an isolation move. Goal of the isolation move is to deform 
the boundary between the reference and the target party 
so that the support of O is sufficiently separated from 
the reference party. See FIG [2] 





T 




() 1 

R 




R 


T 




R 




R 


T 



FIG. 2: Applying the isolation move, the conditional entan- 
glement spectrum is deformed in such a way that i) for the 
new conditional entanglement spectrum, O is sufficiently far 
away from the reference party ii) the difference is a condi- 
tional entanglement spectrum with small conditional mutual 
information. 

Once the support of O is isolated from the reference 
party, we can apply a separation move, which separates 
the support of O from the target parties. See FIG|3] 

Last step is to apply an absorption move. Absorption 
move enables us to write the correction terms as a linear 
combination of Hj^..(j.\g. such that i) the support of O is 



FIG. 3: Applying the separation move, the conditional en- 
tanglement spectrum is deformed in such a way that i) for 
the new conditional entanglement spectrum, O is sufficiently 
far away from both the reference and target parties ii) the 
difference is a conditional entanglement spectrum with small 
conditional mutual information. 



3^ 

R 




R 


T 



R 




H 


T 



cp 

R 




R 


T 



FIG. 4: Applying the absorption move, the conditional en- 
tanglement spectrum is expressed in terms of a linear combi- 
nation of conditional entanglement spectrum HA^.c^\Bi such 
that i) the support of O is contained in either AiBi or Bid 
ii) I{Ai : Ci\Bi) is small. 



contained in either AiBi or BiCi ii) I{Ai : Ci\Bi) — o(l). 
See FIGgl 

To summarize, given a local operator O, one can de- 
compose the conditional mutual spectrum Ha:c\b into 
Ha':C'\b' ^nd correction terms with the following prop- 
erties. First, the distance between A'B'C and the sup- 
port of O is 0{l). Second, the correction term consists of 
sum of conditional mutual spectrum such that the sup- 
port of O is contained in the reference party and one of 
the target parties. Third, the conditional mutual spec- 
tra in the correction term have small conditional mutual 
information for the ground state of topologically ordered 
system. 

In Section |IV| and |V) we shall frequently encounter 
terms of the following form. 



TT{pABcHAi:C,\B,0), 



(21) 



where O is an operator whose support is contained in 
AiBi. If PABC is Co-bounded, this term can be bounded 
as follows. 

Tt(pabcHa-c,\b,0) = TrAiBi^rCi{pABcHA,:Ci\BiO) 

< \^^'Ci{PABcHAi:Ci\B,)\l\\0\\ 

<coI{A:C\B)\\0\\. (22) 
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C. Lieb-Robinson bound 

Lieb-Robinson bound states that there is a constant 
speed of hght so that the correlation decays exponentiaUy 
outside the eflFective hghtcone. We refer the readers to 
Ref E71 for pedagogical introduction to the subject. Here 
we assume the quantum many-body hamiltonian satis- 
fies the Lieb-Robinson bound and study its consequences. 
Given an observable Oa{Ob) supported on A{B), Lieb- 
Robinson bound can be formally stated as follows. 

||[OaW,Ob]|| <c||0^||||OB||min(|A|,|B|)e^i(''*-^(^^^», 

(23) 

where < c,ci,v < oo are some constants that depend 
on the parameter of the hamiltonian and d{A, B) is a 
distance between A and B. 0{t) = e~*^*Oe'^* is a time 
evolution of operator O under the hamiltonian. 

Similar, albeit weaker locality bound holds for 
J'^oo f{^)^A{t)dt if f{t) decays sufficiently fast. To state 
this more formally, we introduce a superoperator defined 
as follows. 

/oo 
e-'"'Oe'"'f{t)dt (24) 
-OO 

It is worth noting that in the energy eigenbasis, 

<ffiO)Uj^ fiE,-Ej)0,j, (25) 

where /(w) is an inverse Fourier transform of f{t). 

We also define a truncated superoperator 4>J by intro- 
ducing a cutoff T. 



Therefore, $^3(7;^) can be approximated by a local op- 
erator. 

Corollary 3. 

\\<i>f,{v^) - [$^^^(^'^)]..Mll < c'\\v^\\e-^^^, (31) 



where Vi{r) is a set of sites whose distance from the 
support of Vi is less or equal to r. 

It would be remiss if we do not mention was origi- 
nally introduced by Hastings in the context of QBP algo- 
rithm and finite temp eratu re correlation decay properties 
of a fermionic system.'^^'^oi j-^ fact, it is not a coincidence 
that these operators appear in seemingly different con- 
texts. As one can see in the following lemma, is a 
quantum channel that appears naturally when comput- 
ing a directional derivative of a density matrix.^ 



Lemma 3. For p{s) 



-liH(s) 



^p(s)|.=o = ^($/f(V^)Ps + /i.c.)-/?($^,(y)), (32) 

where h.c. is hermitian conjugate. 

Similar technique was used by Bachman et al.^^ in the 
studies of the ground state of gapped systems. They 
showed that a unitary evolution that adiabatically con- 
nects the ground states of two different hamiltonian is 
generated by a path dependent generator of the follow- 
ing form. 



(26) 



Lieb-Robinson type locality bound for can be es- 
tablished as follows. 

||[$/(Oa),Os]|| < ||[$J(Oa),Ob]|| + ||[AJ(0^),Ob]||, 

(27) 

where Aj = $y — $J. The first term can be bounded by where T = minggjo^i] r(s) and Writ) is some superpoly- 



D{s) 



D{s)U{s),U{0) = I. 



(33) 



(34) 



dH{s) _ 



V. 



\f{t)\dt\\OA\\\\OB\\ce 



c'(vT-d(A,B)) 



(28) 



from the Lieb-Robinson bound. Second term can be 
bounded by 



-T.T] 



\f{t)\dt\\OA\\\\OB\ 



(29) 



nomially decaying function. In our setting, 
Each of the local terms Vi in V can be approximated as 
follows.l2ll 

W'^wAv^) - [■i>wA^i)]v^(r)\\ < C\\v,\\G^'H^), (35) 

whe re v is the Lieb-Robinson velocity appearing in Equa- 
and G^(a;) is a function that satisfies the foUow- 
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Depending on the function /, one can optimize the 
bound with a judicious choice of T. An example that 
illustrates this idea is /f (w) = *'^"|ji^2^^^ • 

Lemma 2. // H satisfies Lieb-Robinson bound, 

mf40A),0B]\\ < c\\OA\\\\OB\\nim{\A\,\B\)e-^^^^ , 

(30) 

for some constant < c, c' < 00. 



tion 



mg property. 



K 



130e^a;^"M2/7(x) 



< X < a;o 

X > xq. (36) 



Estimates for the constants are K sa 14708, 36057 < 
Xq < 36058.^^ Also, Ua{x) is defined as follows. 



Ua{x) 



(37) 
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IV. GROUND STATE OF EXACTLY SOLVABLE 
MODELS 



Exact formula for the entangleraent entr opy is k nown 
for quantum double and Levin- Wen modelsPEHIII] jf the 
subsystem is simply connected, the entanglement entropy 
satisfies area law. 



Sa = a\dA\ - 7, 



(38) 



where 7 is the topological entanglement entropy. These 
systems have zero correlation length, so the density ma- 
trices of two nonoverlapping regions factorize, i.e. pab = 
Pa® Pb- Therefore, the following formula holds for the 
entanglement entropy. 



S 



AB 



(39) 



if A n B = 0. 

Using standard perturbation theory, for a family of 
quantum states p(s) that are differentiable with respect 
to s. 



dSA rr ^dp~ 



(40) 



Therefore, 

dI{A : C\B) 
ds 



= Tr{^^HA..c\B) 

= lJ2'^Y{[<^Wr{Vj),P^]HA:C\B), (41) 

j 



where Pq is a projector onto the ground state. 

Without loss of generality, let us consider terms Vj that 
are distance al or less away from ABC, where a > 
is some constant. Using deformation moves, Ha:c\b — 
Ha':C'\b- + Y.^a^HA,.c.\B^, whcrc d{Vj,A'B'C') = 0[l) 
and I{Ai : Ci\Bi) = 0. By Petz's theorem, HAi:Ci\Bi = 0. 
Now approximate $i^p(wj) by [^WT{vj)]vj(ci) for some 
c > such that the support of ['^WT{vj)\vj(ci) does not 
overlap with A'B'C . This implies the following relation. 

i:TmwM)]v,ici),PG]HA':C'\B') = 0. (42) 

To see this, consider an operator O that is supported 
on one of A', B' , C , or £> = {A'B'Cy. 



iTr{[0,Po]HA':C'\B') = ;T7Tr(e'"*Poe-^"'-ffA':C'|B 



dt 



(43) 



where the infinitesimal generator generates a unitary 
transformation supported on {A'B'C'y. Since the entan- 
glement entropy is invariant under local unitary transfor- 
mation, this is 0. The correction terms are of the follow- 
ing form. 



Using Equation 16 and [35j we conclude that the effect 
of each terms are bounded by cJG'^^\c'^)Pd for some 
constant c and c'. 

Terms that are distance al or more away from 
ABC can be bounded by approximating ^Wr{''^j) 
[^Wr{'Vj)]vj{R)i where R is the distance between Vj and 
ABC. Combining all of these contributions together, we 
arrive at the following bound. 

%=o<cJ{-ri\yric'-) (45) 
as V V 

for a sufficiently large subsystem size One can see that 
the bound diverges for gapless systems. 

We note in passing that the same technique can be 
applied to topologically trivial configurations, i.e. I{A : 
C\B) = 0. Under general perturbation that consists of 
sum of short-range bounded- norm terms, conditionally 
independent configurations become approximately con- 
ditionally independent. One may wish to establish a 
bootstrapping argument that recursively uses the approx- 
imate conditional independence of these configurations. 
Main difficulty of this approach lies on proving the cq- 
boundedness. 



V. STABILIZER MODELS AT FINITE 
TEMPERATURE 

Unlike the ground state of the exactly solvable mod- 
els, exact formula for the entanglement entropy of a fi- 
nite temperature system is not known except for few spe- 
cial cases.^i^'SSlMI Xo cope with this difficulty, we make a 
nontrivial but natural assumption: that the corrections 
from the deformation moves consist of conditional mutual 
spectrum with small conditional mutual information. For 
3D toric code, topological entanglement entropy does not 
depend on the size of the subsystem for a sufficiently large 
subsystem. tSI We shall denote the conditional mutual in- 
formation in the correction terms as e{l) and study how 
the first order perturbation effect depends on it!^ We 
shall also assume that the correlation decays exponen- 
tially. 



C{Oa,Ob) < ||OA||||OB||min(|A|,|i3|)e 



-d{A.B)/i 



(46) 



Stabilizer model refers to a hamiltonian of the following 
form 



H 



(47) 



(44) 



where Ji > are coupling constants and s^s are elements 
of the stabilizer group. Stabilizer group is an abelian 
subgroup of Pauli group. ^'^ Important examples include 
toric code, color code, their higher dimensional general- 
izatio ns, and other glassy topologically ordered systems 
I|i4 | 35 jt38] j.^pQ^.^Q^j^|- property of the stabilizer mod- 
els is that their reduced density matrices commute with 
each other. 
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J2s es{A) '^i^i /or some coefficients 



Lemma 4. pA 

{C^}■ 

Proof, p can be expanded as a sum of stabilizer group 
elements. After taking the partial trace, any operator 
that has nontrivial support on A'^ vanishes. Any stabi- 
lizer group element that has nontrivial support only on 
A survives. But these terms are generated from the gen- 
erator of the stabilizer group, so they are again elements 
of the stabilizer group. □ 

ft trivially follows that for the Gibbs state of the sta- 
bilizer hamiltonian, reduced density matrices commute 
with each other. Therefore, any reduced density matrix 
PABC for the stabilizer model is 1-bounded. To see this, 
note the following inequality 



i:ii(lnDi - InDa) > Di - D2 



(48) 



for positive semidefinite operators Di , D2 which com- 
mute with each other. Setting Di = pabc £^nd D2 = 
PabPb^ Pbc and taking a partial trace over C, we con- 
clude that Ttc{pabcHa:c\b) is a positive semidefinite 
operator. Since li norm is equal to the trace for positive 
semidefinite operator, pabc is 1-bounded. 

Consider terms Vj that are distance al or less away 
from ABC. Using the deformation moves, Ha:C\b = 

Ha':C'\B' + T.^a^HAr.c.\B^, wherc d{v„A'B'C') = 0{l) 
and I{Ai : Ci\Bi) = e{l). Choose an approxima- 
tion radius R such that ^jii {vj ) is approximated by 

['^ ff>{vj)\vj{R)- First order effect of Vj can be divided 
into three parts: the connected correlation between 
[$^3 and Ha':C'\B', the approximation error of 
^ ff and the corrections from the deformation moves. 
Terms that are distance al or more away from ABC can 
be similarly bounded by using the exponential correla- 
tion decay and making a judicious choice for the approx- 
imation radius R. All of these effects combined together 
results in the following bound. 



=0 < 0(Z'^(e-^i'/«)+0(?2^e-^^'/^))+0(/^e(/)), 

(49) 

where D is the number of spatial dimensions and ci , C2 , C3 
are some numerical constants. 



VI. COMMENT ON HIGHER ORDER TERMS 

Close inspection of the first order bound reveals that 
the CQ-boundedness plays a pivotal role in the derivation. 
For example, consider a perturbed ground state of the 
topologically ordered system which satisfies area law ap- 
proximately. Equation |45] is only modified by including 
the area law correction terms, provided co-boundedness 
is guaranteed. 

ft turns out that the CQ-boundedness in a finite neigh- 
borhood of s implies a nontrivial bound for the higher or- 



be applied to topologically trivial configuration as well 
as topologically nontrivial configuration. Since Equation 
l45l relied on the fact that the conditional mutual infor- 
mation of topologically trivial configuration is small, we 
can bootstrap this argument to the bound higher order 
terms. 

Assuming the co-boundedness for s £ [0, sq), following 
inequality holds. 



ds 



I{A : C\B),\ < Ss{l) + : C,|B,) 



(50) 



where Ss{l) is a function that decreases superpolynomi- 
ally with I, and is a finite number that is uniformly 
bounded for s S [0, sq]. I{Ai : Ci\Bi)s is a conditional 
mutual information appearing in the correction terms of 
the deformation moves. 

ft the energy gap remains open for s e [0, sq), ^s(0 
can be uniformly bounded by some S{1) that decays su- 
perpolynomially in L As a result, one can obtain the 
following recursive bound. 



l7.-7o|< / 6il)+y2adiA,:C^\B,)s'ds' 
Jo 

J Jo Jo 

(51) 



der terms as well. The key idea is that Equation 45 can 



Here we used the fact that the conditional mutual infor- 
mation arising from the deformation move is at s = 0. 
Recursively applying this logic, the second order term 
can be bounded by 0{P6{1)). Higher order terms can be 
obtained in a similar manner. 

To investigate the validity of the CQ-boundedness for 
general quantum many-body system, we have generated 
random density matrices and studied a relationship be- 
tween both sides of Equation [18] The result is plotted in 
FIGj5j For pure states, we have applied a random uni- 
tary from Haar measure. For mixed states, we have ran- 
domly generated eigenvalues from uniform distribution 
over [0, 1], normalized, and applied random unitary from 
Haar measure, ft seems that for certain states that have 
small conditional mutual information, smallest value of 
Co increases significantly. For this reason, we urge the 
readers to be careful in using this condition in general. 
This difficulty can be circumvented for stabilizer mod- 
els against stabilizer perturbations, since commutativity 
of the reduced density matrices is preserved. However, it 
remains to be seen if the correction terms from the higher 
order deformation moves are small. 



VII. CONCLUSION 

We have demonstrated that conditional independence 
strongly constrains the structure of quantum many-body 
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I(A : C\B) 

FIG. 5: We have numerically computed I{A : C\B) and 
\T^rcpABcHA:C\B\i for 10*' randomly generated pure states. 
Largest observed ratio \Ttc PabcHa:C\b\i/ I{A : C\B) was 
24.21924. 

system so as to ensure the first order perturbative sta- 
bility of the topological entanglement entropy. Admit- 
tedly, our technique gives bounds in limited settings 
where i) exact conditional independence is achieved or 
ii) reduced density matrices commute with each other. 
However, once these conditions are met, the argument 
can be applied quite generally. In particular, we ex- 
pect our method to be applicab le to t he studies of Cha- 
mon's model and Haah's model.'^^'^ These models sat- 
isfy the topological quantum order conditions introduced 
by Bravyi et al, and their hamiltonian consists of sum 
of frustration-free commuting projectors.^! Therefore, the 
energy gap is protected against generic perturbation that 
consists of sum of short-range bounded-norm terms. 

There are compelling reasons to believe that these 
models are not described by BF theory or multiple stacks 
of Chern-Simons theory: movement of the quasiparticles 
are constrained in a peculiar manner, and their ground 
state degeneracy is determined by some number-theoretic 



function that depends on the size of the system.lSMo] j^. 
would be interesting if one can apply our method to find 
a linear combination of entanglement entropy that allows 
the first order perturbative stability. 

We have also shown that our method can be ex- 
tended to higher orders of perturbation series if the cq- 
boundedness holds in a finite neighborhood of s, but such 
statement seems unlikely to hold for general quantum 
states. It would be very interesting if one can find an 
alternative technique that relies on the conjecture intro- 
duced in Ref|28l There author was able to show that 
the connected correlation function between conditional 
mutual spectrum and local operator vanishes if the local 
operator is supported on one of the subsystems, provided 
certain extension of strong subadditivity is true for gen- 
eral quantum states. Unfortunately, local operators that 
are supported on multiple subsystems are bound to ap- 
pear, as shown in the analysis of the present paper. 

As for the finite temperature topological entanglement 
entropy in 3D, we needed two nontrivial assumptions to 
bound the first order perturbation effect. First, the con- 
nected correlation function between two observables de- 
cay exponentially. Second, the correction terms from the 
deformation moves can be expressed as a sum of small 
conditional mutual information. We emphasize that nei- 
ther of these assumptions were explicitly proved. Fur- 
ther studies in explicitly bounding both of these terms 
are necessary. 

While the structure of conditionally independent state 
is relatively well understood, much less is known about 
the states that are approximately conditionally indepen- 
dent. We hope our work motivates further studies in such 
direction. 
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